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EFFECT OF TEMPERATURE ON THE ELECTRONIC INSTABILITY 
AND THE CRYSTALLINE PHASE CHANGE AT LOW 
TEMPERATURE OF V 3 SI TYPE COMPOUNDS 


J. Labbe and J. Frledel 

College of Sciences, Solid State Physics, 91 Orsay 


Introduction 


In a preceding article [ 6 ], we showed that a Jahn-Teller /303 

effect on the d band structure can explain the instability of the 
cubic phase of type V 3 SI Intermetallic superconductor compounds 
at zero temperature, with the tetragonal phase being stable. In 
this second part, we study the effect of temperature. By applying 
Fermi statistics to our model, we show that the cubic phase re- 
covers its stability above a certain temperature. This allows 
the martensite transition observed at low temperature by Batterman 
and Barrett [4] to be. understood. In addition, we show that, in 
our model, the transition very likely is of the first order but, 
meanwhile, it is' only accompanied by slight discontinuities in 
the variation of different physical parameters, such as lattice 
distortion and the elastic constants. Likewise, the latent heat 
of transformation certainly is very low. 

In the second section, we stress more particularly the vari- 
ation of the elastic constants with temperature, and we show that 
the experimental results obtained by Testardl, Bateman, Reed and 
Chlrba [3] really appear to reveal a first order transition. 


In the third section, we apply our calculations numerically 
to' the example of 'V 3 SI. 

1. Effect of temperature and order of the transition 


For a finite value of the temperature T, we have to find the 
value of e which gives the minimum free energy F=U-TS. The Fermi 
level E'p is a function of e and T. Its value is given by the 
equation 
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where Ej?(0), Is the Fermi level for e=0 and T=0. At a nonzero 
temperature, no band Is ever completely empty and we cannot dis- 
tinguish the cases e<ec j and-eJ^'E- 

The free energy is a function of c and T; for a given temper- 
ature T, its variation with e is written: 
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where Ep is the Fermi level at temperature T and e=0, and- recall- 
ing that there are 3 N transition atoms in all in the crystal 
and Q electrons per transition atom in the d band under consider- 
ation. 'In (II. 2), the variation 'of the contributions to entropy 
with e, other than that of the d electrons, is disregarded. 

Equation (II. 1) can be solved numerically by machine . cal- 
culation. By successive approximations, E'p is calculated for 
various values of e and T. The corresponding numerical values 
of dF are then obtained from (II. 2). In this way, a system of 
curves is obtained, which respresents, variations of dF with e 
for various temperatures. Fig. 1 shows the system obtained in 
this way for ■V 3 SI, as we shall see in Section 3. It is seen that 
the cubic phase (e= 0 ) again becomes stable above a certain tem- 
perature Tjjj. • 

In conclusion, at very low temperatures, the electrons re- 
main confined in the bottom of the band, where the density of 
states is very high, which involves instability of the cubic phase 
On the other hand, when the temperature rises, the band filling 
is modified, to the benefit of a lower density of states, which 
reduces the effect of electron instability. The presence of the 
term 1/2 A'e^, which changes very little with temperature, thus 
permits the cubic phase to regain its stability. To determine 
the order of the transition, we have to find the detailed shape 
of the curves which correspond to neighboring temperatures of Tm, 
for very small deformations e. Then, dF can be expanded with 
respect to e, and the order of the transition can be found, in 
such a way that the coefficients of this expansion vary with tem- 
perature. Therefore, if the following is written: 
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Fig. 1. Variations dF of free 
energy with distortion e for 
various temperatures T. 

Key: a. In eV per vanadium atom 


coefficients A, B and C, ob- 
tained from (II.l) and (II. 2) 
by calculations too long. to be 
reproduced here, are expressed 
by: 
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and 


f(E) is the Fermi . function 
f" (E) , . . . are its successive derivatives with respect 
E. y is the coefficient of the term in e^ in the 

level, and it is 


X iS 
displacement 


E'p-Ep of the Fermi 


f'(E), f"(E), 
to energy 
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B' and C' are the anharmonic contributions of the conduction elec- /305 
trons and nontransition atoms; B' can be roughly estimated from 
the equation of state of Griineisen. He finds B’=^-6A' (B'<0). It 
is very difficult to evaluate C, but we shall see that it is 
reasonable to disregard it in a certain temperature region. 
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A machine calculation permits A 
various temperatures T. Figs. 2k, 
behavior of their variations with T, 
V3Si. 




Fig. 2. Variations with tem- 
perature T of coefficients A, 

B and C of the expansion 

dF = 1 1^.4 

of the free energy for very small 
distortions . 

Key: a. In eV per vanadium atom 
b. T in. °K 


, B and C to be determined for 

2B and 2C also indicate the 
for the particular case of 

It follows from the 
theory of Landau [1] that 
the transition can only be 
second order, if coefficients 
A and B cancel for the same 
temperature. In principle, 
this can occur for a well 
determined value of the 
initial filling of the band. 
But, this possibility can 
only be considered accidental. 
It would be quite extraordin- 
ary for the filling to have 
just the required value. On 
the other hand, group theory 
[2] shows that B has no cause 
to be identically zero due 
to crystal symmetry. There- 
fore, if the change in the 
structure under study here 
does not make any other , 
parameter than uniform dis- 
tortion e Interfere the 
transition is certainly of 
the first order. But, we 
find that, in fact, in the 
cases of interest to us., 
coefficients A and B cancel 
for similar, temperatures, 
so that, if the transition 
really is of the first order, 
it is, meanwhile, only ac- 
companied by slight discon- 
tinuities in variation of 
the physical parameters. In 
addition, it only brings a 
low latent heat into play. 


We call To the temperature for which- A cancels by changing 
sign. Generally, B and C are zero for this temperature, and the 
transition is first order. In the immediate vicinity of Tg, A 
can be expanded in T-Tg, and A-K(T-Tg) can be written, where 

K is a positive constant. If this vicinity is small enough, B and 
C can be identical to their values at Tg, or Bo=B(Tg) and Co=C(Tg). 
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Cq has to be positive for there to be a transition. The number 
of solutions of the equation dP(e)=0 depends on the sign of the 
discriminant 


A = IbI-~^CoK{T-To) 

which cancels out for the temperature 

T,., = To +{2li)){nUC,K) (II. 7) 

and there then is a double root 

c(7’„) = --(2/3) (Z/o/Co). (II. 8) 


The first order transition occurs at temperature Tm. Bq can 
be positive or negative, according to the initial filling of the 
bands. Figs.- 3A and 3B show the behavior of the variations of dF 
with e for very' close temperatures Tin, respectively, in the cases 
Bq> 0 and Bq<0. • 



.The- latent heat of trans- 
formation is equal to, the, vari- 
ation AU of the internal energy 
between the states {T=Tm, e=0} 
and {T=Tm, e=e(Tm)). The . 
Helmoltz relation permits the 
writing of 

At; = AF— 7’(a/D7’) (AF). (II. 9) 

The free energy is the same 
for the two states under consid- 
eration. Therefore ,. AF=0 . On 
the other hand, the following is 
obtained from (II. 3) 


Fig. 3- Variations dF of free 
energy with distortion e for 
temperatures T very near trans- 
ition temperature Tjj], in .the 
cases Bq> 0 (Fig. 3a) and Bq<0 
(Fig. 3b) . ■ 
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■- .j-Si-C^V-srAK 


But, since dF has a minimum at e-e(Tni) > we have •‘. 

A{T„) c{7’™) + B{T„) t\T„) + C(T„) = 0 

and, therefore, 

c*{T„) . 

(II. 10) 

' 2. Variation of elastic constants with temperature 

In our model, we can predict the behavior of the variations 
with temperature of the total modulus of shear Ap which is con- 
nected to the constants of rigidity C 1 1 and C 12 , calculated for a 
crystal volume containing 3N transition atoms, by 

3Abl, = ^ (C„ — C,,). ■ • (II. 11) 

For T>Tm., dF is at a minimum for e = 0, and Ap is the coef- 
ficient A, given by equation (II. ^). Contribution A' of the 
conduction -electrons is positive, and it varies very little with 
temperature. On the other hand, contribution A" of the d electrons 
expressed by' the first two terms in (II.^), is negative, and it 
decreases very rapidly in absolute value as the temperature in- 
creases. For temperatures immediately above Tm, A’ and A" have 
almost equal absolute values, and Ap is very small. On the con- 
trary, at ordinary -temperatures , A" has an absolute value of no 
more than one tenth of the value of A'. Ap is then nine tenths 
of A’ . Fig. 2A shows that, just above Tjn, the variations of Ap 
with T can be considered linear, and their extension to the left 
Intersects the temperature axis at T=Tg, with Tq very slightly 
less than Tj^ (Fig. ^). Difference Tj^-Tq Is obtained from (II. 7). 
This anomaly in the behavior of Ap has been brought out experi- 
mentally in 'V' 3 Si by Testardl, Bateman, Reed and Chlrba [3]. 

For T^Tm, dF is at a minimum at a certain value e(T) of e, 
different from zero. Then, Ap is no longer equal to coefficient 
A given by (II. 4), but it. is directly connected to the second 
order coefficient of the expansion of dF in e-e(T). Thus, it is 
found that 

At = A -\-2Ut{T) +3CzliT) + fJ)z^{T) (11.12) 

where D is the coefficient of the fifth order term in (II . 3) , 
which can be written De^/5. 

For T=Tm, (II. 8) and (II. .12) show that At (Tm) =A(Tm) + ^De ^ (T^), 
which means that the modulus of shear has adlscontinulty for DE^(Tm’) 
during the transition though it is very small, because it is pro- 
portional to E^(Tm). This is quite conslstant with the fact that,' 
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Tm 
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the transition Is of the first 
order , 



Pig, 4. Variations of modulus 
of shear At with temperature 
T. ■ ■ . 


3. Application to case of V 3 SI 


In V3SI, the distortion 
observed by Batterman and Barrett 
[4] is tetragonal, with a positive 
value of e. Therefore, we are 
in b or c of the discussion of 
I,2.E [6], and the number Q of 
electrons in the nearly empty d 
band per transition atom satisfies 
the double inequality 




<<?< 




which, with (1,11) expanded by 
Q, is also written 


Morin and Malta [5] have measured the electron specific /307 

heat ln'V3Si, at temperatures between 0®K and 2p°K, and the value 
they obtain for the density of states at the,, Fermi level, can be 
considered as relating to the tetragonal phase. They find a 
numerical value of 11 eV“* per vanadium atom, with the two direc- 
tions of spin taken into , account , The following then is obtained 
from (1,45) 

= 3A'X11- 


or, by disregarding the small term'aq0/2. 


i^Fi) X 11 ^ 25 cV-f 

On the ether hand, Batterman and Barrett find 

a(l +e)A.(l-|) = 1,0025 

or e=l . 67 • 10" ^ , 


The following then is obtained from (1,11) and (1,21) 

1,07 X 10-V~ 8«y/--yr«(/iK)- (II. l4) 

Experimentally ,• at ordinary temperatures,' it is found that 
[3] Cii= 2. 87*10^^ erg‘cm-^ and C 12 = 1 . 20 • 10 erg*'cm- ^ . 

The crystal parameter is 2a=4.72 A. The number of vanadium 
atoms per cm^ is, therefore 3N-6(2a)"^ =5. 72*10^^. 

From •( II . 21 ) , it is then found that, at 300°K, At = 27.4 eV per 
vanadium atom. According to Section 2, it is then seen that 
there has to be A'=;30 eV per transition atom. 


7 






Relation (II. 1^) involves aq="1.57. This leads to acknowledge- 
ment, for the unknown parameter, which is the coefficient of 
Slater, of the vanadium atom in the crystal, the quite reasonable 
numerical value q = 0.67 A~^'. 

The width 2|,Efj,| of the d band under consideration is, like- 
wise, a parameter which we no not know, but, from (11.13), it is 
found that 

1,82 cV < < 2/,2 eV. 


This is a very reasonable order of magnitude for the nar- 
rowest d band, i.e., the d 2 ^^x 2 _y 2 band. 

For example, let us take \2 Eui|=2.20 eV. Prom (1. 11), Q=0.0593 
is then found. It is found in case b of r.2.E. 

The machine calculations reported in 1 then give the following 
results, which can be, read in Figs. 2 A,B and C: 

K = 1 eV degree"^ per vanadium atom 
Bq=-0.5*103 eV per vanadium atom • 

Co= 1.6*10® Ev per vanadium atom 
Tq = 13.3°K. 

(II. 7) and (II. 8) then give 

.Tm=13.4°K and e (T^) =0 . 2 • 10- ^ 

The experimental value of Tj„ is clearly greater than 20 to 
25°K, Actually, while our unidirectional model really results in 
a high, narrow peak for the density of states, it does not give 
us its exact shape in the immediate vicinity of the bottom of the 
band with certainty. But, for a very slightly occupied band (here, 
Ep-Em-lO-® eV), this shape begins to be of Importance, as soon as 
the temperature rises a dozen degrees above absolute zero. 

Fig. 5 represents the calculated variation of e with temper- 
ature T. It is seen that, when T Increases, initially, e hardly 
varies. On the other hand, when T approaches Tjn, e begins to de- 
crease very rap.ldly . This appears to be in agreement with experi- 
ment [4]. . 

For T=Tui, e only has an eighth of'lts value at 0°K. There- 
fore, the transition is only ‘ accompanied by a slight discontinuity, 
which is quite difficult to detect experimentally. Furthermore, 
the calculation gives (Fig. 2 ) : 
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per vanadium atom, 

ft tr !t 


It t! 


It 


Formula (II. 10) then gives 
ation AU-0.006 calorie per gram 
value is difficult to measure. 



Fig. 5. Calculated variations of • 
distortion c with temperature T 
in the tetragonal phase. 

Key: a. T in °K 

be expected. The absence of a peak 
[ 5 ] is in agreem.ent with a first or' 


for the latent heat of transform- 
atom in transition. Such a low 


Finally, since the trans- 
ition is first order, hyster- 
esis can be expected. With 
the Inhomogeneities of the 
crystal taken into account, 
this undoubtedly spreads out 
the transition temperature. 
Different portions of the 
crystal do not undergo all 
the changes of structure at 
the same temperature. Only 
measurements made on a speci- 
men which is very suitable, 
from the metallurgical point 
of- view have a chance of . 
making an infinite specific 
heat peak appear, which cor- 
responds to the existence of 
latent heat. On the other 
hand, if the transformation 
was second order, a specific 
heat peak of finite height, 
but without hysteresis, would 
n the measured specific heat 
r transformation. 


To summarize, among the parameters which interfere in our 
calculations, crystal parameter a and the coefficient of shear A' 
at ordinary temperatures are known experimentally with certainty. 
On the other hand, the. specific heat measurements only give us an 
indication, of the order of magnitude of the density of states n(Ep) 
at the Fermi level. Likewise, in reference to the values usually 
excepted for the transition elements, we have only a rough idea 
of the width 2|Em| of the narrowest d band and of the coefficient-, 
of Slater q of the vanadium atom in the crystal. 

By using these 5 numerical values, our calculations permit 
Interpretation of the low temperature values and the thermal 
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variations of two physical quantities: distortion in- the 
tetragonal phase and a coefficient of shear. Distortion e at 
low temperature has a value, which is in good agreement with that 
measured by Batterman and Barrett. They also allow the difference 
in electron specific heat in the two phases, tetragonal and cubic, 
to be understood. ' ’ 

Transition temperature T^ is only obtained in order of mag- 
nitude. Actually, different numerical tests show that it only 
varies slowly with band width. However, it depends greatly on 
details of the density of states n(E), in the vicinity of the 
Fermi level, and these are little known.' 
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